Anisotropic Mean-Field Theory for Semiflexible Polymer Dynamics under Tension 
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We introduce an anisotropic mean-field approach for the dynamics of semiflexible polymers un- 
der tension. The theory is designed to exactly reproduce the lowest order equilibrium averages of a 
stretched polymer, and includes the non-trivial influence of long-range hydrodynamic coupling. Vali- 
dated by Brownian hydrodynamics simulations, the theory is highly accurate over a broad parameter 
range, even far from the weakly bending limit. It captures the change in viscoelastic response due 
to polymer stiffening under tension that has been observed experimentally in cytoskeletal networks. 



Understanding semiflexible biopolymers under ten- 
sion is relevant both in vivo — determiningthe mechan- 
ical response of cytoskeletal networks [H H, 9 — an d * rl 
vitro, due to a proliferation of single-molecule techniques 
like optical tweezers that involve manipulating stretched 
DNA 0. Yet for one of the most basic questions there 
is still no comprehensive, quantitatively accurate theory: 
what is the dynamical response of a semiflexible chain 
prestretched by an external force? Existing approaches 
are typically based on the weakly bending approximation 
( WBA) [5( , which assumes the polymer contour remains 
nearly straight — strictly valid only in the limit where the 
persistence length l p is much greater than the total length 
L, or for a large force F ^S> ksT/lp. This has proven a 
highly useful simplification, widely applied in both equi- 
librium [f| HI, 0] and non-equilibrium [t| contexts 
for stretched semiflexible polymers. Yet such theories are 
of limited value in describing the complex crossovers be- 
tween dynamical regimes occurring for weaker forces or 
more flexible chains. Moreover, most of these theories 
ignore long-range hydrodynamic interactions, or include 
them only at the level of assuming distinct longitudinal 
and transverse friction coefficients, £ii w Cl/2, appro- 
priate for a rigid rod 0. This only renormalizes times 
scales, without affecting the dynamic scaling. As we will 
argue below, long-range interactions do have a significant 
influence on the dynamics even in the weakly bending 
limit, and must be incorporated to obtain a quantitative 
comparison with experiments. 

In this work, we address the gaps in the existing ap- 
proaches by constructing an anisotropic mean-field the- 
ory (MFT) for a semiflexible chain under tension, in- 
cluding hydrodynamics through a pre-averaging approx- 
imation. The solvent-mediated coupling between all 
points on the polymer is not an incidental element in 
the method, but a key to its success: the long-range 
interactions make the mean-field approach more realis- 
tic. In fact, at F = 0, isotropic MFT ll] is able to 
capture precisely the dynamics of an end-monomer in 
DNA strands observed through fluorescence correlation 
spectroscopy [12j . Without any fitting parameters, it 
gives an accurate description for strands 10 2 — 10 bp 
in length, covering a wide range between stiff and flex- 
ible chains [HI]. Unlike earlier (isotropic) attempts to 
extend MFT to F ^ 0, EE El] , our anisotropic the- 
ory is designed to yield the exact lowest-order equilib- 



rium properties for directions parallel and perpendicular 
to the applied force, which are derived from the numerical 
quantum solution of the worm-like chain (WLC) model. 
Having fixed the correct static quantities, we show that 
the theory gives highly accurate and non-trivial predic- 
tions for the dynamics, validated through comparisons 
with Brownian hydrodynamics (BD) simulations. In par- 
ticular, it reproduces the change in viscoelastic proper- 
ties due to semiflexible polymer stiffening under tension, 
which has been observed experimentally in cytoskeletal 
networks put under stress either through deformation jy] , 
or the activity of motor proteins [![. 

We begin by reviewing the isotropic MFT approach 
to stretched semiflexible polymers [lj, EE Ell , and high- 
lighting its limitations. The starting point is the WLC 
Hamiltonian, {7wlc = (lpksT/2) J ds(d s u(s)) 2 — Fz ■ 
J dsu(s), which describes the elastic energy of a space 
curve r(s) with contour coordinate < s < L and tan- 
gent vector u(s) = d s r(s) constrained by local inexten- 
sibility to |u(s)| = 1 Vs. The first term is the bending 
energy, parametrized by the persistence length l p , while 
the second term is the external field due to a force F along 
the z axis. The |u(s)| = 1 constraint makes the dynamics 
of the system analytically intractable, but the partition 
function Z, expressed as a path integral over u(s), can 
be approximated through the stationary phase approach, 
yielding a Gaussian mean-field model. The end result 
is Z ss cxp(— /3-7-mf) = / fu exp(— (3Umf), where lo- 
cal inextensibility has been relaxed and the MFT Hamil- 
tonian is: Umf — (e/2) / ds (d s u(s)) 2 + v j dsu 2 (s) + 
Uq (u 2 (0) + u 2 (L)) — Fz ■ j ds u(s). From the stationary 
phase condition, O^J-'mf — 9v (j ^mf = 0, one finds that 
the parameters v and vq are functions of F and e, and 
act as Lagrange multipliers enforcing the global and end- 
point constraints J ds(u 2 (s)) = L, (u 2 (0)) = (u 2 (Z>)) = 
1. For F = Q, Umf Las been studied extensively (l7l.fla|. 
and setting e = (3/2)/ p fcsT it reproduces exactly the 
WLC tangent-tangent correlation (u(s) • u(s')) and re- 
lated quantities. For F ^ 0, it successfully yields the 
known asymptotic forms for the average end-to-end ex- 
tension parallel to the force: (R z )/L = 2l p F /iksT for 
F -> 0, and (R z )/L 
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^W^r/^TpF for F -> oo 

[Til ], where R = ds u(s). These agree with the Marko- 
Siggia exact result for the WLC [l!| , except for the factor 
3/8, which should be 1/4. 

Underlying this discrepancy is a serious problem in the 
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FIG. 1: (Color online) For a polymer with L/l p — 5, the end- 
to-end vector fluctuation (5R 2 ) and its components, (SR 2 ), 
(SR 2 ± ), with varying F (derived from the exact quantum solu- 
tion of the WLC), compared to the isotropic MFT prediction. 
Inset: cloud of end-point positions (red points) taken from a 
BD simulation of a polymer (N = 50 beads), with L/l v = 5 
and force F — 20 ksT/lp applied along the vertical axis. The 
initial polymer configuration is shown in blue. 



isotropic MFT: for large F it cannot correctly account for 
the anisotropic fluctuation behavior of the WLC. To illus- 
trate this, let us consider displacements 6R\\ = R z — (R z ) 
and SR± = R x or R y . The MFT does not differentiate 
between || and _L fluctuations: (SR. 2 ) = (SR 2 ^) for all 
F. However for the WLC, (SR.?,) becomes much smaller 
than (SR 2 ^), as is evident in the inset of Fig.[TJ showing a 
snapshot of end-point fluctuations taken from a BD sim- 
ulation. A similar anisotropy exists between the _L and 
^fluctuations at F — for a stiff filament with L < l p 
[20l | . Using the exact mapping of the WLC onto quan- 
tum motion over the surface of a unit sphere [2l[ , one can 
numerically calculate the fluctuation magnitudes, which 
are shown in Fig. [TJ (For details see Sec. 1 of the sup- 
plementary material [22|.) Since the MFT averages over 
all coordinate directions, its estimate for the magnitude 
converges to the exact (SR 2 ^) at large F, as the _L fluctu- 
ations dominate in that regime. It misses entirely the || 
component. Thus to understand the dynamic response of 
stretched semiflexible polymers one needs a more suitable 
theoretical starting point. 

To resolve these difficulties, we propose an anisotropic 
version of the Gaussian model: 



U = E y / ds (d s u a (s)f + v a J ds u 2 a (s) 
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Wc now have a total of 



seven parameters: {e a , v a , voa}, ol =||)-L, and a force 
rescaling factor \. The guiding philosophy will be sim- 
ilar to the isotropic case: to construct a Hamiltonian 
which closely approximates the equilibrium behavior of 
the WLC under tension, and then use it to derive a dy- 
namical theory. For this purpose, we require that U 



should reproduce exactly the following lowest-order WLC 
averages: (R {{ ), (SR 2 a ), J Q L ds (u 2 a (s)) , (4(0) + u\ (L)>. 
The latter can all be calculated from the quantum solu- 
tion to the WLC [22| , while the corresponding analytical 
expressions derived from U lead to seven equations for 
the seven unknown parameters. These can be solved nu- 
merically for any given L, l p , and F (examples are shown 
in Sec. II of 22]). In the limit F — ► 0, our approach recov- 
ers the stationary phase condition of the F = isotropic 
model, as expected. 

Our dynamical theory builds on the hydrodynamic 
pre-averaging approach used earlier for the isotropic 
MFT [ll|, [23[. Here we give a brief outline of the ap- 
proach, with the full details in Sec. Ill of '22] . The time 
evolution of the chain r(s, t) is governed by the Langevin 
equation: r a (s,t) — — J ds' £i" vg (s — s')SU / Sr a (s' ,t) + 
£ a (s,t), where £ a (s,t) are stochastic velocities, and hy- 
drodynamic effects are included through the pre-averaged 
anisotropic mobility /i"vg( s — s ')- The- latter is derived 
from the continuum version of the Rotne-Prager ten- 
sor ^Jjt(s, s';x) [Hj], describing solvent-mediated interac- 
tions between two points s, s' on the contour at spatial 
separation x. If the equilibrium probability of finding 
such a configuration is G(s, s';x), then the integration 
J d 3 x'jt(s, s'; x)G(s, s';x) yields a diagonal 3x3 tensor 
whose a =||,-L components we denote as /i" vg (s — s'). 
In the absence of hydrodynamic effects (a case we will 
consider as a comparison), the free-draining mobility is 
Hfd(s — s') — 2afioS(s — s'), where a is a microscopic 
length scale (i.e. the monomer radius), and /iq is the 
Stokes mobility of a sphere of radius a. We assume 
the stochastic velocities £(s,i) are Gaussian, with cor- 
relations given by the fluctuation-dissipation theorem: 

(Z a (s Metis',?)) = 2k B TS(t - i')Mavg(s - «')■ _ The 
Langevin equation, together with boundary conditions 
at the end-points due to the applied force, can be solved 
through normal mode decomposition, yielding all the dy- 
namical quantities which we will analyze below. 

To validate the anisotropic MFT, we compared the 
theoretical results to BD simulations of a bead-spring 
worm-like chain, consisting of N spheres of radius a (con- 
tour length L — 2aN) , with the monomers hydrodynam- 
ically coupled through the Rotne-Prager tensor. (Full 
simulation details are in Sec. IV of 22]). Fig. O shows 
results for a representative semiflexible polymer, with 
L = 100a and L/l p = 5. For each direction a, the mean 
squared displacement (MSD) of the end-to-end vector, 
Af(t) = ((R a (t)-R a (0)) 2 ) [Fig. (Ha)], and an end-point 
of the chain, A c a nd (t) = ((r a (L,t)-r a (L,0)) 2 ) [Fig.^b)], 
is depicted at two different forces F. There is excellent 
quantitative agreement with the BD simulations (dashed 
curves), with the maximum errors w 10% for the _L and 
w 20% for the || results in the time ranges shown. The 
biggest discrepancies occur at short times for the || com- 
ponent with F = 20 ksT/lp, where the length scale of the 
motion is comparable to the bead size, and we expect the 
discrete BD chain to deviate from continuum MFT be- 
havior. The close agreement is all the more remarkable 
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FIG. 2: (Color online) For a polymer with L = 100a, L/l p = 
5: (a) A£(t); ( b ) A™ d (t); (c) ImJ^H; (d) ImJ™ d (w). Solid 
lines are the anisotropic MFT results, while dashed lines are 
taken from BD simulations. In all cases results are given for 
a =_L, || at two forces, F — 2 and 20 ksT/lp. Filled circles 
mark the relaxation times r a i, derived from the MFT, while 
the insets show the local slopes of the MFT curves in the 
log-log plots. 
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3: (Color online) (a) 
100a, L/l p — 5, F -- 



A" ld (t), a =_L, ||, for a chain with 
= 20 knT/lp. The top two curves 



include hydrodynamic interactions, while the bottom two are 
free-draining, (b) A® e (£) for a chain with L = 100a, L/l p = 
1/3, F = 60 k B T/l p . In both (a) and (b) MFT results are 
shown as solid lines, BD simulations as dashed lines. The red 
curve in (b) is from the combined MFT+WBA theory. The 
insets show the local slopes of the MFT curves in the log-log 
plots. MFT relaxation times T a i are marked by circles, while 
the crossover times t* are marked by vertical dashed lines. 



since the MSD shows a complex crossover behavior that is 
only imperfectly captured by asymptotic scaling theories, 
for example the Granek prediction [5| for _L motion based 
on the WBA. For t -C Tli, the longest relaxation time of 
the chain in the _L direction, Granek derives two regimes 
separated by the crossover time t* = 2l p ksT '/3F 2 hqcl: a 
stiffness-dominated regime at t <C t*, with MSD cx i 3 / 4 , 
and a force-dominated regime at tj_i 3> t 3> t* , with a 
slower scaling cx t 1 / 2 . The insets of Fig. HJa)-(b) show 
the local slopes of the log-log MSD plots, dlog A Q /dlogi, 
with times tj_i calculated from the MFT marked by dots. 
With increasing F, we do indeed find the local slope is 
reduced, but the dynamic scaling is modified by two im- 
portant effects: (i) the slow crossover to center-of-mass 
motion at times t t±i, where the slopes of A° nd and 
A^ e approach 1 and respectively; (ii) logarithmic cor- 
rections due to hydrodynamics, which increase the local 
exponent on the order of 10%. Note that even in the 
strongly stretched limit, Flp/ksT 3> 1, where the poly- 
mer is nearly straight, hydrodynamic effects are signifi- 
cant. Fig.[3Ja) shows MFT and BD results for A% nd with 
and without hydrodynamics for a chain where L = 100a, 
L/l p = 5, and F = 20 k B T/l p . The MSD components 
in the two cases cannot be related through a simple time 
rescaling: for f < t < r a i, we see clearly the expected 



i 1 / 2 behavior for the free-draining chain (local slopes are 
shown in the inset), while the exponent is pushed up to 
Rs 0.6 — 0.7 with hydrodynamics. 

Through the fluctuation-dissipation theorem, the MSD 
can be used to extract the viscoelastic properties of 
the chain: Fourier transforming the MSD functions 
gives the imaginary parts of the end-to-end and self re- 
sponse functions of the polymer end-points, Im J£(lu) = 
—iuA^(uj)/2kBT, [i = ee,end, which are defined as 
J™ d M = Sr a (L,u)/f a (u), J£(w) = SR a (uj)/f a (uj). 
Here 8r a {L,uj) and 8R a (u>) are the complex oscilla- 
tion amplitudes resulting from a small force f a (io) — 
/oexp(— iuot) applied to one or between both ends of the 
chain in addition to the prestretching tension F. Using 
the MFT, one can express J£(oj) as a sum over normal 

mode contributions, J£(w) = J2 n =i ^n/f 1 ~ wr an ) + 
Sp^cndiDa/uksT, where the center-of-mass diffusion pa- 
rameters D a , relaxation times T an , and coefficients A£ n 
are determined by the MFT solution [22j|. The mode 
number cutoff M = L/8a is chosen to roughly model 
the discrete nature of the chain at length scales compa- 
rable to the bead diameter, but results at larger length 
scales are independent of the cutoff Tl]. Fig. [DJc)-(d) 
shows Im J° c (oj) and Im J™ d (w), for the same parame- 
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ters as in Fig.[5Ja)-(b), compared to the results extracted 
from BD simulations (additional plots of the real com- 
ponents are in Sec. V of [13] )■ The good quantitative 
agreement with BD in the time domain is carried over 
to frequency space: the simulation trends are accurately 
reproduced by the MFT. For lu -C , we see a mainly 
elastic end-to-end response, J° c w ^ai(l + iu>T a {), with 
an effective spring constant (A®®) -1 . The self response 
of the end-point at these small frequencies is propor- 
tional to the center-of-mass mobility, J™ d « iDjujk B T. 
For lu > we pass into the more interesting high- 
frequency regime governed by the complex nature of nor- 
mal mode relaxation under tension and hydrodynamic 
interactions (up to the ultraviolet cutoff at t~ m , above 
which the discreteness of the chain dominates). The ef- 
fects of tension in this regime have been directly observed 
in cytoskeletal networks through microrheology p], Q: 
with increasing force the dynamic compliance is reduced, 
and the high-frequency scaling changes from lj~ 3 / 4 (the 
behavior of a relaxed semiflexible network) to lu^ 1 / 2 . 
We find both of these stiffening effects in our MFT re- 
sults, with an important correction: unlike a network, 
where hydrodynamics is screened, in the single polymer 
case the long-range interactions modify the local slopes: 
rather than —3/4 and —1/2, we see w —0.8 at weak 
force changing to —0.6 at strong force. This correc- 
tion, along with the full crossover behavior of the imag- 
inary response — proportional to the power spectral den- 
sity (PSD) — should be observable in future nanorheology 
experiments for single semiflexible chains (i.e. the AFM 
techniques already used to extract the PSD of flexible 
polymers US!)- 

Though we have demonstrated the accuracy of the 
MFT, one must be careful to consider the limitations 
of our approach. For small forces, Fl p /ksT <C 1, and 
extremely long chains, L 3> 1%/or, one expects self- 
avoidance to be important [26| , which is entirely left out 



of the MFT theory. (For DNA with l p w 50 nm, ex- 
cluded volume effects would appear at L » f 2 ^m.) In 
the opposite limit of extremely stiff chains, with l p ^> L, 
we find another problem: for a system that is nearly a 
rigid rod, no Gaussian model will be able to capture the 
|j dynamics. In Fig.^b) we show A° e (i) for a chain with 
L = 100a, L/l p = 1/3, and F = 60 k B T/l p . We see that 
the MFT underestimates the || relaxation time by an or- 
der of magnitude, saturating to equilibrium much quicker 
than the BD result. However the _L Gaussian theory is 
still remarkably precise, deviating < 7% from the BD 
curve throughout the entire time range. This provides 
us with an alternative solution: use the _L MFT results 
in combination with WBA to estimate the || quantities. 
The key relation is Un(s,t) w 1 — Uj_(s, i), valid in the 
stiff limit. Using the _L MFT estimate of Uj_(s,i), we 
have derived a second-order perturbation expansion for 
Ajj e (£) (details are in Sec. VI of [H), yielding the red 
curve in Fig. [3Jb) . This gives a much better agreement 
with BD (< 25% deviation). Thus with the help of WBA, 
we can extend the scope of our Gaussian MFT even to 
the rigid rod limit. 

In summary, we have developed an anisotropic MFT 
for the dynamics of semiflexible chains under tension, 
whose most notable features are an unprecedented quan- 
titative accuracy — verified through BD simulations — 
and a broad range of validity, not restricted to the 
weakly bending regime. Understanding kinetics of single 
stretched chains is interesting in itself (or as the first step 
toward more elaborate theories of stressed networks), but 
it can also be exploited in other contexts: optical tweezer 
force-clamp experiments depend sensitively on the dy- 
namical response of the DNA handles that are attached 
to the object of interest, whether a nucleic acid hairpin or 
protein [27j • The implications of our theory in analyzing 
such experiments will be pursued in future work. 
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Supplementary Material for "Anisotropic Mean-Field Theory for Semiflexible 

Polymer Dynamics under Tension" 

Michael Hinczewski and Roland R. Netz 

I. QUANTUM SOLUTION OF THE WLC 

The mapping between the WLC and a quantum particle moving on the surface of a unit sphere can be exploited to 
calculate exactly many of the equilibrium properties of the system [l[ . Here we follow a technique similar to Ref. Q 
to numerically evaluate thermodynamic averages of the WLC to arbitary accuracy. To compute all the quantities of 
interest, we start with a WLC Hamiltonian augmented by two extra terms: 

u= lpkBT r dg ( dsU ( s)) 2 _ F f dsUz ^_ F •[ dsu x (s)~K [ dsu 2 z (s). (S.l) 



The extra parameters F x and K will later be set to zero after taking the appropriate derivatives of the partition 
function. Rescaling the contour variable as r = s/l p and factoring out fcgT, we can rewrite U in a simpler form: 



/3U = dr 



T 







\{d T n{r)f - fu z (r) - f x u x (r) - ku\(r) 



(S.2) 



where T = L/l p , f = f3l p F, f x = (3l p F x , and k = j3l p K. Let us define the propagator G(uq, uj; T) as the path integral 
over all configurations with initial tangent u(0) = u and final tangent u(T) = u^: 



G(u ,u T ;T)=/ PuJJ^(u 2 (s)-l) exp (-/?[/) . (S.3) 



Then for boundary conditions with free end-point tangents the partition function is given by Z = 
(47r)~ 2 J s duo duT G(uq, ut; T), where the integrations are over the unit sphere S. The quantum Hamiltonian corre- 
sponding to (3U is 

H = -(1/2)V 2 - /cos^- / x sin6»cos0- fccos 2 <9, (S.4) 

describing a particle on the surface of a unit sphere. In terms of the associated quantum eigenvalues E n and eigenstates 
xp n (u), the propagator G is given by: 

G(u , u T ; T) = J2 e- E " T r n (uo)MuT) = £ e-^o^a^l^UoJlWur) , (S.5) 

where in the second step we have expanded out the eigenstates in the basis of spherical harmonics, ip n (u) = 
Ylim a n£m^m(u), with coefficients a n i m . For a given n, these coefficients are just the components of the nth 
eigenvector for the Hamiltonian H in the Y\ m basis. Thus to proceed, one needs the detailed form of this ma- 
trix: Hi,rn;l>,m' = W?, m;I , , m / + n{ m .y m , + W^.,, + H k lm . Vm ,. We list below only the nonzero elements of each 
contribution that are relevant to the calculation, with symmetry of the matrix implicitly assumed: 

W? )0i j,o = Kl + 1). Ko;i+i,o = -f( l + !)[( 2 ' + !)( 2 ' + 3 )] _ 1/2 > 
/k _ k(2l(l + !)-!) ^ fe fc(/ + !)(/ + 2) 



n l. 0:1,0 ~ A,n , i\ o - ' 'N.Oi 



4Z(i + l)-3 ' w+2 '° (2Z + 3)V4Z 2 + 12Z + 5' (S.6) 



_1. / (/ + m + l)(Z + m + 2) , _1 / (f-m+l)(i-m + 2) 

■: | + 1 . m + 1 ~ 2 /x V 4Z(Z + 2) + 3 ' >Wn,m-i- 2 /*y 4 /(/ + 2) + 3 



The eigenvectors a n z m and eigenvalues E n can be readily calculated numerically by truncating the infinite matrix 
Hi,m;i',m' to a finite size (with cutoff chosen large enough to get the desired precision). The partition function Z can 
then be written as: 

Z = T4-T2 / du ° duT e ~ E " Ta ni>m' a ni m Y l f m ,(u Q )Y lm (u T ) = — ^ e~ EnT a* nQO a nOQ (S.7) 

* ' S n,l,m,l',m' ^ n 



Most of the thermodynamic averages used in the anisotropic MFT can be directly derived from Z: 
d 

-AogZ 

/ so =0,fc=0 



(i?ll) = l P f- f \ogZ 



ds{u\{s))=L- [ ds{ul(s))= Ip-^logZ 
o Jo a & 



/ x =0,fc=0 



J /x=0,fc=0 

The end-point averages are calculated similarly: 



(5Rl) = 2ll^\ogZ 



/ x =0,fc=0 



(ug„ (0) + «g„ (£)) = 2 - (^ ± (Q) + ?i 2 ±(L)) = _JL_ jf duQ duT U 2 |G(UQ) Ut . t) = £ e-iVTa 



(S. 



a n00 , a n2Q 



67T 

(S.9) 



II. CALCULATING THE PARAMETERS OF THE ANISOTROPIC MFT 



With free end-point tangent boundary conditions, the partition function 

Z=(Att)- 2 I du du L / X>uexp(-/3J7) 

JS -/u(0)=u o 



(S.10) 



corresponding to the anisotropic MFT Hamiltonian U [Eq. (1) in the main text] can be evaluated analytically: 
2 v / 2^ 3/2 /?- 1 e^ ± 



Z 



(e> 2 + 4j/ 2j sinh ( LcJ± ) + 4e_Li4u.au cosh(Luu) \ ^a w a + 4z/ 2^ sinh ( Lw|| ) + 4ep ||W|| cosh(LW||) 



^IKi 



• exp 



> X 2 ,F 2 (sinh (^) ( £ ||Lc 2 - 4^„) + 2Lv n ^ cosh 



(S.ll) 



2e||w| f e||W|| sinh f + 2i> || cosh 



where w Q = ^J2v a /e ai a =||,_L. Similarly one can extract analytic expressions for all seven of the thermodynamic 
averages used in the fitting of the MFT parameters: 



(#11 



/ ds(wf(s)) = -^- 1 ^-log^ / dsK 
Jo °^ll Jo 



(9 2 



(SRl 



Llu i 



o 

4^ j 



2i/ ± coth (^) + e ± u ± J ' 



( M 2 || (0) + M 2 || (L))=-/3- 1 



d"o\\ 



logZ, 



(v»±(P)+i%j_{L)) = -l3- 1 



\ogZ. 



(S.12) 



By setting these expressions equal to the exact WLC results calculated from the quantum approach, Eqs. (|S.8[) and 
(|S.9j) . one obtains a system of seven coupled equations that can be solved numerically for a given L, l p , and F, yielding 
the seven Hamiltonian parameters: en, ej_, vu, v±, Vqu, vo±, and \- Fig. IS - 1 1 shows a set of solutions for L = 100a, 
l p = 20a, and varying F. In the small force regime, F <C ksT/lp, where stretching is negligible, the parameters 
converge to the same values as in the isotropic MFT: z/|| = u±_ = iksT / Al p , i/ || — u o± = 3fcsT/4, e|| = e±_ = 3kBTl p /2 
II- 



III. DYNAMICAL THEORY OF THE ANISOTROPIC MFT 



Here we adapt the dynamical theory used successfully for the isotropic MFT [_| [_| to the anisotropic case of a chain 
under tension. The general hydrodynamic preaveraging approach is similar to that used for the Zimm model 0, [H • 
The time evolution of the chain r(s, t) follows the Langevin equation: 



/ ds' */x* (s, s'; r(s, t) — r(s', t)) 
Jo 



SU 



5r(s',t) 



+ €(«.*)■ 



(S.13) 
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FIG. S.l: Parameters of the anisotropic MFT Hamiltonian U as a function of force F for a chain with L — 100a, l p = 20a, 
derived as described in Sec. I and II in order to reproduce exact equilibrium averages of the WLC. To plot all parameters in 
one graph, y-axis units have been chosen such that ksT = a = 1. 



Here the £(s, t) is the stochastic contribution, and *Ji(s, s'; x) is the continuum version of the Rotne-Prager tensor [f|, 



*fi(s,s';x) = 2a^5(s - s') 1 + Q(x - 2a) 



8tti]x 



X (g) X 



x~ 



4iri]x J 



(S.14) 



describing long-range hydrodynamic interactions between two points at s and s' on the chain contour, separated by 
a spatial distance x. The microscopic length scale a corresponds to the monomer radius, r\ is the viscosity of water, 
/io = 1/QTTija is the Stokes mobility of a sphere of radius a, and the function excludes unphysical configurations 
(overlap between monomers). 

Eq. (|S.13|) cannot be solved directly, because the hydrodynamic tensor depends on the exact configuration of the 
chain at time t, so we employ the pre-averaging approximation: replacing < Ji(s,s';r(s,t) — r(s',t)) with an average 
over all equilibrium configurations, /i avg (s, s'): 



Vavg(S)S') = J d 3 x^t(s,s';x)G(s,s';x) 



(S.15) 



where G(s,s';x) is the equilibrium probability of finding two points at s and s' with spatial separation x. For the 
anisotropic Hamiltonian U this probability takes the form: 



G( S , S ';x) 



2ira±(s — s') \2irau(s — s') 



1/2 



exp 



3xi 



2a±(s-s') 



3(X|| -F|s-s'|/2^||) ; 
2<7n(s-s') 



(S.16) 



where a a {l) = (3(|i|cJ Q + exp(— — l)//3e Q w^,. In deriving G we have assumed a large chain length i, which 

simplifies the resulting analytical expression. Plugging Eq. (|S.16[) into Eq. (|S.15[) leads to: 



^ vg ( s - o o o 

MavgO," 8 ') = I Mavg( s - s/ ) 

A*L g (s - s')j 



(S.17) 
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where the anisotropic mobilities /x" vg can be written in terms of integrals over coordinates X — ■ ^1 x ^ I X jj and £ = 

3 3 / 2 e(Z-2a>o /- 1 ^-3ttC 2 , 3 ,x,, lW \ (C 2 - l) -x 2 3(p£ — Fl/2v\\) 2 \ 
Mav g (0 = 2a/j 5(i) 

3 3 / 2 6(^ - 2a) Mo f 00 , f 1 , a tt (-3C 2 (x 2 - 2) + 9a; 2 - 2) f 3 (( 2 - l) x 2 3((x - Fl/2v«) 2 \ 

1 ax I dQ — 1 exp 1 1 



(271)3/2^(0^0 Aa J-i 4x ^ 2a ± (l) 2a l{ (l) J' 

(S.18) 



These integrals are evaluated numerically to obtain the mobilities as a function of contour distance I. In Fig. IS. 21 
we show the results for L = 100a, l p — 20a, and F = 1.0 kgT/a. Note that the mobility parallel to the stretching 
direction is enhanced relative to the transverse component, as we expect for an extended chain. 
The pre-averaged version of the Langevin equation can now be written as: 

£r a (M) = - [ ds'^ vs (s- S ') f +Z a ( s ,t), (S.19) 
at J s 5r a {s',t) 

for a =||, _L. The £(s,i) are Gaussian random vectors, whose components have correlations given by the fluctuation- 
dissipation theorem: 

(Us, *)£«(*', 0) = 2k B T6(t - t')^ vg (s - s'). (S.20) 

Plugging in the form of U , the internal force term in the Langevin equation becomes 

SU d 4 d 2 

r ( i i\ = ^—M 8 '^) - ^ a —^r a (s',t) ee 0%r a (s',t), (S.21) 

OV(x\S j tj US US 

where we have introduced the differential operator Of, . To complete the dynamical theory, we must specify the 
boundary conditions at the chain ends: 



d 3 ,„ , „ d ,„ N d 3 d 

; -—rr a {L,t) +2v a — 
as 6 as 



" ea ^3 r «( '*) + 2^ a — r Q (0,<) = FS aj \\, -e a -^r a (L,t) + 2v a —r a (L,t) = FS a ^ 



d 2 d d 2 d (s ' 22) 

e a Tr^r a (0,t) - 2v Oa —r a {0,t) = 0, -e a - -~r a {L, t) - 2v Qa —r a {L,t) = 0. 

OS z OS OS* OS 

The first two represent the force applied at the ends, while the second two the absence of torque. To properly deal 
with the boundary conditions for F ^ 0, we write r(s,t) in the following way: r a (s,t) = r a (s,t) + F6 a ^(j)(s), where 
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f (s, t) satisfies the homogeneous (F = 0) version of Eq. (|S.22[) . while <p(s) is chosen such that the total function r(s,t) 
satisfies the full boundary requirements. The resulting form for <p(s) is: 



Ls Lisas' 2 Vo\\ 



2(^11+2^11) 2e|| (Lz/y + 2^ || ) 3e||(Lz/|| +2z/ ||)' 



(S.23) 



We now proceed to transform the Langevin equation into matrix form, which will allow us to solve it through 
numerical diagonalization. Let us assume £ a (s,i) satisfies similar boundary conditions to f Q (s,f), and expand both 
functions in normal modes ip"{s), with amplitudes p an {t) and q a n(t) respectively: 



00 00 



r a (s,t) = 2^Pan(*)<(s), 6*(M) = L«»(^"( S )- ( S ' 24 ) 



n=0 n=0 



The normal modes ipn(s) are chosen to be eigenfunctions of the differential operator O", satisfying 0"ip"(s) 
^an*Pn( s ) f° r eigenvalues X an . These eigenfunctions take the form [f|: 



,,, , C an ( T r sinK an s sinhG Q „s . ,c.>-, 

^ {S) = V— r°" co8g an L/2 +Ga " coshG a „L/2 j ' n ° dd ' (S ' 25) 



;a/ x C an ( r , cosies cosh G an s \ 

= V~ l~^" sin^„£/2 +G "" sinhG a „L/ 2 J ' 

where 

G^ n — if^ n = 2v a / e a , Aq = , A Q n = ta-R'an + 2j/ q G 2 n . (S.26) 

The eigenfunctions obey the boundary conditions in the F = version of Eq. (|S.22[) , which fixes the constants K an 
and Gam while the C an are normalization coefficients. Plugging the normal mode expansions into the Langevin 
equation, and exploiting the orthonormality of the Eqs. (|S.19[) - (|S.20[) become: 

dt ^0 ( S - 27 ) 

(q an (t)q am (t')) = 2k B T5(t - t')H% m , 

where 

H nrn = \ dsj ds' ^° (s)/Cg( S ~ S ')^Z( S ') > 

(S.28) 



L 



ds 



ds' V ) |(s)MLg( s ~ sO 2 ^!? 11 L , o ■ 

6||(il/|| +2l/ ||) 



Both -ff " m and tu n can be evaluated through numerical integration. In order to make solving these equations feasible, 
we introduce a high-frequency cutoff M on the mode number, keeping only the slowest-relaxing modes n — 0, . . . , M— 1, 
whose hydrodynamic interactions are described by the leading MxM sub-blocks of the matrices H a . Following Ref. Q 
we set M = L/8a, which provides good agreement at short times with Brownian dynamics simulations of bead-spring 
chains with monomer radius a. For longer times, where the polymer motion is on length scales much larger than a, 
the dynamical results are insensitive to the precise value of the cutoff. 

The final step in simplifying the dynamical theory is diagonalization. Let J a be the MxM matrix with elements 
Jn m = H" m X am , A Q „ be the eigenvalues of J a , and C a the matrix diagonalizing J a : [C a J a (G a ) _1 ] nm = A an 5 nm . 
Assuming non-degenerate eigenvalues A arl , the matrix C a also diagonalizes H a through the congruent transformation 
[C a H a {C a ) T ] nm = QanSnm, defining parameters 6 Q „ [7[. The diagonal version of Eq. (|S.27[) then reads: 

d 

Ql P ^n{t) = -A an Pan(t) + FW n S a ,\\ + Qan(t) , , g ^ 

(Qan(t)Qam(t')) = 2k B TS(t - t')5 m ,n®an, 
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where 



M-1 



M-i 



M-1 



Pan(t) — J]] Cn m Pam(t), Qan(t) — C" m q am (t), W n — CjJ 



(S.30) 



and r a {s,t) = E„^m(*)*£0) + ^a,||0(«) with modified normal modes *£(s) = £ m <(s)[(C t, )" 1 ]mn- Using 
Eq. (|S.29[) it now becomes possible to solve for a variety of dynamical observables. For example, the result for the 
MSD of a chain end-point is: 



AT(t) = {(r a (L,t) - r Q (L,0)) z ) = 2k B T 



e 



n>0 



2D a t + 2fc s rJ] A™ d (l - exp(-A Q „t)), 



(S.31) 



n>0 



where we have introduced the center-of-mass diffusion constant D a = ksTQao^^ (L)) 2 , and coefficients A e ^ 
8 a „(f"(L)) 2 /A Qrl . Similarly, for the MSD of the end-to-end vector, 



A~(<) = ({Rait) - R a (0)) 2 ) = 2k B Tj2 " exp(-A Q „t))«(L) - ^(0)) 



n>0 



An 



2k B T Cnt 1 ~ expHW)), 



(S.32) 



n>0 



where A- = 6 Qn (*«(i) - 4<£(0)) 2 /A Q 



IV. BROWNIAN DYNAMICS SIMULATIONS 



For the BD simulations [sj used to validate the mean-field theory, the chain consists of N beads of radius a whose 
positions Ti(t) are governed by the discrete Langevin equation: 



d,Yi(t) ( 8U(r 1 ,...,r N ) 

2^ M v 



dt 



3 = 1 



(S.33) 



Long-range hydrodynamic interactions between monomers are included through the Rotne-Prager [lfj| mobility matrix 
/x ij : 



A 4 



=A*o<5j. 



1 



(S.34) 



where r ^ = r j — . This matrix also determines correlations for the Gaussian stochastic velocities ^ i (t) according to 
the fluctuation-dissipation theorem: 



(€i(t)®€,-(0) =2k B T*j?ijS(t-t'). 



(S.35) 



The elastic potential of the chain U = Uben + C4tr + Ulj + U ey ± consists of four parts: (i) a bending energy 
t^bcn = (e/2a) 5Zj(l — cos^i), where 9i is the angle between two adjacent bonds, and e is related to the persistence 
length l p as e = l p k B T; (ii) a harmonic stretching term f/ str = (7/4a) (r i+1 j — 2a) 2 where inextensibility is 
enforced through a large modulus 7 = 2000fcaT/a; (iii) a truncated Lennard- Jones interaction J7l,t = w Si<j ©(2a — 
rij)[(2a/rij) 12 — 2(2a/rij) 6 + 1] with ui — 3k B T; (iv) an external force F along the z direction, J7 ex t = —Fz ■ (rjv — ri). 

In the numerical implementation of Eq. (|S.33|) . the Langevin time step is r = 3 x 10~ 4 a 2 /(k B T[i ), where /j,q is 
the Stokes mobility of a monomer, and a typical simulation lasts ~ 10 8 — 10 9 steps. Data is collected every 10 2 — 10 3 
steps, and averages for the dynamical quantities discussed in the main text are based on 5-25 independent runs. 
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V. RESULTS FOR REAL AND IMAGINARY PARTS OF RESPONSE FUNCTIONS 



1 1 1 L I — I — I — 1 1 I I I I I 1 1 




w [k B Tno/ll] 



FIG. S.3: For a polymer with L = 100a, L/l p = 5: (a) Re J c a c {to); (b) Re JT d {to); (c) Im J^w); ( d ) Im J^H- Solid lines 
are the anisotropic MFT results, while dashed lines are taken from BD simulations. In all cases results are given for a =_L, | 
at two forces, F = 2 and 20 ksT/lp. Filled circles mark the relaxation times r Q i, derived from the MFT, while the insets show 
the local slopes of the MFT curves in the log-log plots. 



VI. WEAKLY BENDING APPROXIMATION FOR Af{t) BASED ON THE 1 MFT RESULTS 

For small deviations from the rigid rod limit, the || and _L tangent vectors of the WLC can be related as: u« (s, t) — 

yl — Uj_(s, t) w 1 — Uj_(s,t)/2 — Uj_(s,t)/8 + • • • , where = {u Xl u y ). This is the fundamental equation for the 
weakly bending approximation, and it allows one to derive certain aspects of the |j dynamics assuming the _L dynamics 
are known, specifically the behavior of Uj_(s, t). As seen in Fig. 3(b) of the main text, the anisotropic MFT provides 
a highly accurate prediction for the _L end-to-end MSD even for very stiff chains, so we can exploit the reliability of 
the _L dynamical theory through the WBA approach. 

We focus on finding an estimate for the || end-to-end MSD function Ajj c (t), though the method is generalizable 
to other dynamical quantities. Ajj c (£) can be expressed as Ajj e (£) = 2(C||(0) — C^(t)), where the correlation function 
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Cm (t) is given by: 



q|(t) = (i?n(t)i?n(0))-(ii||) 2 = / & / d s '("n(s,t)w||( s ',o)) 



L 

<is (s, 0)) 



(S.36) 



Here we have used the fact that R\\(t) = J Q dsu^(s,t). Plugging in the second-order expansion M||(s, t) ~ 1 — 
u 2 L (s,t)/2 — Uj_(s,t)/8, we get an expression for C\\(t) involving averages over various products of u 2 L (s,t). From the 
normal mode expansion in Sec. Ill we know that u±(s,t) — d s r±(s,t) — J2 n P±n(i)^n'( s ); where Pj_„ = (P xn , Pyn) 
and ^„'(s) = 9 s ^/„(s). Thus all averages over Uj_(s,i) are averages over the normal mode amplitudes Pj_ n (t), and 
these can be directly calculated from Wick's theorem and the solution of Eq. (|S.29|) for the _L components. The final 
result for C»(t) at order 0(u.j_) has the form: 



I 

C \\( f ) =Z)/fc(*)/l(*) M W + 4 E fk(t)fl(t)f m (0)M kl M mmkl + ]T / fc (W)/m(*)/n(*)M, 



fe.J k.l.m k,l.\ 



where: 



4 fk(t)fl(t)fm(0)fn(0)M kl 

rain, an • 



fk{t) = — t exp(-A_L fc t), 



(S.37) 



M kl = [ ds^i'(s)^'(s), (S.38) 
Jo 



Jo 

Thus C||(t) can be expressed entirely in terms of quantities from the _L MFT solution: the parameters {Aj_„, 0_i_ n } 
and the normal modes {'F^(s)}. Numerical evaluation of Eqs. (|S.37[) - (|S.38j) yields C\\{t) and hence Ajj e (i). 
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